Neutrino oscillations: Measuring #13 including its sign 
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In neutrino phenomenology, terms in the oscillation probabilities linear in sin #13 lead naturally to 
the question "How can one measure #13 including its sign?" Here we demonstrate analytically and 
with a simulation of neutrino data that Pep and P MM at L/E — 271-/A21 exhibit significant linear 
dependence on 613 in the limit of vacuum oscillations. Measurements at this particular value of L/E 
can thus determine not only #13 but also its sign, if CP violation is small. 
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INTRODUCTION 

To incorporate neutrino oscillations the standard 
model is conventionally extended by adding a mass term 
and a mixing matrix. This theory of three- flavor neutrino 
oscillations has been successful in accommodating the re- 
sults of neutrino oscillation experiments, save LSND 0. 
The theory contains six independent parameters: three 
mixing angles 6jk, two independent mass-squared dif- 
ferences Ajk ■= m 2 j — m\, and one Dirac CP phase 6. 
Ref. 0, for example, provides a summary of the current 
knowledge of the values of these parameters and describes 
future experiments. The long baseline (LBL) K2K ex- 
periment is in the process of confirming the results 
of the Super-K atmospheric experiment |^ by measur- 
ing the parameters A 32 and 823, independent of atmo- 
spheric neutrino flux models. Future LBL experiments, 
MINOS 0, OPERA 0, and ICARUS 0, will improve 
upon the bounds for these parameters. Additionally, a 
global analysis of these LBL experiments could provide 
a lower bound on the magnitude of 6*13 [8(. Future LBL 
experiments might also resolve the question of mass hier- 
archy and the level, if any, of CP violation in the neutrino 
sector 0. 

We here examine the related question of how to best 
measure #13, including its sign. In Ref. ^(J, we have 
shown that, even in the presence of matter effects [llj . 
neutrino oscillations can be uniquely and completely pa- 
rameterized with the following bounds on the angles: the 
CP phase 6 lies in the range [0, 7r); #13 lies in the range 
[— 7r/2, 7t/2]; and the remaining mixing angles lie within 
the first quadrant. This choice of bounds has two ad- 
vantages. First, present experiments limit #13 to a small 
asymmetric region about zero 0, 0|. Other choices 
would break this region into two disconnected regions. 
Secondly, the CP violating phase is restricted to the first 
two quadrants; this range is sufficient to characterize all 
CP violating effects. Terms proportional to cos<5 are thus 
able to uniquely determine its value, assuming knowledge 
of all the other parameters. 

In the next section, we analytically examine the terms 
of the neutrino oscillation probability formulae that are 



first order (linear) in #13. These terms are proportional 
to either sin 5 or cosJ, as indicated in Refs. |2j,|l4|. It has 
been suggested [l4|,ll5j that the presence of such terms, in 
part, can explain the excess of electron-like events in the 
Super-K atmospheric experiment 0. In this work, we 
find that experiments which lie in the oscillatory region 
for the small (solar) mass-squared difference, an L /E on 
the order of 10 4 m/MeV, are sensitive to the linear #13 
terms. We further find that by judicious choice of the 
value of L/E the effects of the CP violating phase S can 
be suppressed if S is near zero or tt, while at the same 
time the effect of the linear term in #13 is maximized. In 
the subsequent section, we utilize a simulation 13] of the 
existing neutrino oscillation data to further examine the 
ability of new data to determine #13, including its sign. 
We assume that either CP is conserved or that our choice 
of the value L/E has provided sufficient suppression of 
the CP violating terms. In the final section we summarize 
our conclusions and provide some thoughts on needed 
future theoretical work. 



FORMAL ANALYSIS 

In this section we provide explicit analytic expressions 
for three neutrino oscillations valid for the incoherent 
limit of the atmospheric mass-squared difference. We 
confine our discussion to vacuum oscillations. For LBL 
experiments through the earth, we indicate which val- 
ues of L and E yield the cleanest measurement of #13 
by avoiding significant contributions from matter effects 
. Additionally, we indicate qualitatively the conse- 
quences of straying outside these energies and baselines. 
We use these analytic expressions to examine where the 
effects of the linear terms can best be seen. As the mag- 
nitude of #i3 and the mass-scale ratio a := | A 2 i |/| A 32 1 
are known to be small, one may expand the oscillation 
probability formulae about these parameters. In these 
perturbations (cf. (l6jl. terms which are linear in $13 are 
suppressed by a factor of a ~ 0.03. From this, one might 
conclude that effects relevant to the sign of 9i 3 are for- 
ever relegated to the realm of the unobservable. Here, 
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this is not the case as we look beyond the valid region of 
these expansions. 

We use the standard representation 01 01 the three- 
neutrino mixing matrix with the notation Cjk = cos9jk, 
Sjk = sinBjk, and 5 is the CP violating phase. In a three- 
neutrino theory, the probability that a neutrino with rel- 
ativistic energy E and flavor a will be detected a distance 
L away as a neutrino of flavor (3 is given by 

P a p(L/E) = 5 a0 

3 

-4 ^e(U aj U* k U pk U^) sin 2 ip jk 

j<k 
j,fc = l 

3 

+2 lra(U a3 U* k Uf3 k U^)sin2ip jk , 

j<k 
j,k = l 

(1) 

where tpj k :— Aj k L/4E with Aj k := to 2 — rrvj.. 

Examining the terms which are linear in sin #13 moti- 
vates us to consider the limit in which the oscillations due 
to the mass-squared differences | A.32 1 = A.31 1 ~ 10~ 3 eV 2 
are incoherent while the oscillations due to IA21I ~ 
10 _5 eV 2 are still relatively coherent. In this limit, we 
may take 

sin 2 (^23 = sin 2 9313 = ^, sin 2(^3 = sin2(pi 3 = 0. (2) 

The oscillation probabilities V efi , "P w , and V^t, in the 
limit of incoherent atmospheric oscillations, are then 
given by 

'Pefj, = [| sin 20i2 cos 20i2 sin 20i 3 ci3 sin 

+ sin 2 20i2C 2 3 (c 2 ,3 - s 2 3 S23)l sin 2 <pi2 

+ isin 2 20i3S^3 + 2Jsin2 V 'i2, (3) 

Vf+fi = 1 - [sin 20i2 cos 20 i2 sin 20i 3 ci 3 sin 20 2 3S2 3 C5 
+2 sin 20i2 cos 20i 2 si3 sin 20 2 3 cos 20 2 3C<5 
+ (1 - sin 2 20 12 c 2 )s 2 3 sin 2 20 23 
+ sin 2 20i2 (c 2 3 - s 2 3 s^ 3 ) 2 ] sin 2 (p 12 
— 2c 13 s 23 (l — c 13 s 2 3), (4) 

V^r = [sin 20i2 cos 20i 2 si 3 (f + s 2 3 ) sin 2023 cos 2023C.5 
(1-sin 2 20i 2 c 2 )s 2 3 sin 2 20 23 
-ism 2 20 12 (l + s 2 3 ) 2 sin 2 20 23 
+ sin 2 20i 2 s 2 3 ] sin 2 y> 12 

+ ism 2 20 2 3ct 3 + 2Jsin2^ 12 , (5) 

with 

J = -| sin 20i2 sin 20i 3 ci3 sin 20 23 Si. (6) 



As we are interested in the sign of 0i3, we isolate those 
terms that are odd with respect to 0i3 

Ve»(0m) - 7V-0i 3 ) = 4J sin 2(^i2+ 

sin 20i2 cos 20i2 sin 2013C13 sin 2023C,5 sin 2 ipi2, (7) 

r^(e 13 ) - r^(-e 13 ) = -4( c 2 3 - S 2 3S 2 3 ) 

x sin 20i 2 cos20i 2 Si3 sin 2023C5 sin 2 (fx 2 , (8) 
?V (#13) - 7V(-0is) = 4 J sin 2<p 12 + 2sin20 12 
x cos20i 2 si3(I + s 2 3 ) sin 2023 cos20 2 3C,5 sin 2 9512. 

(9) 

Note that the sign of 0i3 exhibits the maximal effect 
whenever sin 2 9312 is maximal for c$ ~ 1. This occurs 
whenever cpi2 — (2n + l)7r/2 or, in other terms, for 
L/E = 2(2n + I)7r/Ai2. These oscillations will be more 
coherent for the smaller values of n. 

This choice is fortuitous in that it also removes from 
Eqs. J7J) and terms dependent on s$ as s'm2ip 2 i = 
whenever 9521 is an odd-integer multiple of n/2. This 
removes the CP violating terms from consideration. The 
remaining terms are modulated by c$. Thus for 5 near 
zero or n, we would have a clean measurement of 0i3. 
Also, should CP violation be found to be maximal, then 
the terms involving eg vanish. 

A consequence of removing the dependence on J is 
that CP violating effects are suppressed. At these local 
values of L/E we have V a /3 — P-gn-, where a indicates 
an antineutrino of flavor a; assuming CPT is invariant, 
this can be expressed as V a p = Vp a - This supports our 
previous statement that it is sufficient to only consider 
V efJ , 1 'Pftft > and "Pftx m regards to their dependence on the 
sign of 013. These probabilities in addition to P ee , which 
is a function of 2 3 , will give us all the other oscillation 
probabilities at this value of L/E. The remaining prob- 
abilities are 

Ver = 1 - Vee - P e „, (10) 
V T T = l-V IXT -Ver, (11) 

so that the dependence of V eT on the sign of i3 can be 
surmised from the statements made concerning "P eM and, 
likewise, the behavior of V TT can be surmised from V^ T 
and Vefj,- In what follows, we will assume that any CP 
violation is small so that we may set the phase equal to 
0. 

We have previously performed (l2T | a simulation, as- 
suming no CP violation, of the world's neutrino oscilla- 
tion data. An analysis which includes more recent data 
[l3|. in preparation, produces A% 2 := \ 2 ~ Xmin as pic- 
tured in Fig.^ Included in the analysis are data for neu- 
trinos from the sun [l^, from reactors 0] .atmospheric 
neutrinos [23, and beam-stop neutrinos For one 

standard deviation, the analysis bounds 0i3 to lie within 
[—0.17,0.22] with two minima located at 0i3 = 0.11 and 
—0.04. For the absolute minimum 0i3 = 0.11, we find 
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FIG. 1: The value of A\ 2 = X 2 ~ Xmin versus 613 as ex- 
tracted from the world's neutrino oscillation data using the 
simulation of Ref. |T^I . All other parameters other than #13 
are varied. There are two minima, #13 = 0.11 and —0.04. The 
dashed line represents the results when the "one mass-squared 
dominance" approximation is utilized. The horizontal lines 
represent the 90% and 95% confidence levels. 



12 = 0.48, 23 = 0.80, A 21 = 7.7 x 10~ 5 eV 2 , and 
A 32 = 2.6 x 10~ 3 eV 2 . 

Since Vee is a function of sin 2 #13, the asymmetry seen 
in Fig. [I] must arise from the atmospheric and K2K data, 
which involve V e ^ and V^. If we employ the "one mass- 
squared dominance" approximation, as is often done, we 
find the dashed curve presented in Fig. ^ This approxi- 
mation gives oscillation probabilities that are a function 
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of sin 2 

In order to demonstrate the relative size of the effect 
of the sign of 6*13, we choose some realistic values for the 
mixing angles: 0X2 = 0.56 and 6*23 = 0.78. The first two 
peaks of sin 2 <px2 occur around L/i?=1.6xl0 4 m/MeV 
and 4.8 x 10 4 m/MeV. For such values of L/E, the os- 
cillations due to A 3 2 and A3X would be incoherent. It is 
clear from Eq. © that the screening effect of maximal 
mixing for 023 results in independence of the sign of 0x3 
for V^t- For the remaining two oscillation channels, we 
have a sizable effect. The oscillation probabilities evalu- 
ated at the one-standard-deviation points for 0x3, —0.17 
and 0.22, are 



V, 



-0.17) = 0.35, V eil (0x3 = 0.22) = 0.50, 



(12) 

P w (0i3 = -0-17) = 0.37, 7^(0x3 = 0.22) = 0.22. 

(13) 

The relative differences are more appropriate quantities 
to consider; for we have the most significant effect 
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-0.17) 



-0.25, (14) 
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FIG. 2: Pee versus L/E. The solid curve corresponds to 
the mixing parameters for a fit with #13 = 0. The dashed 
curve represents the two curves given by #13 = ±0.2 which 
are identical. 



while the effect is still large for V e)Jj 

?V(0i3 = 0-22) - 7y(0i3 = -0.17) 
Pe M (0X3 = 0.22) + 7%(0X3 = -0.17) 



0.17. (15) 



We compare the one-sigma extremes of 0x3 in order to 
demonstrate the potential size of the effect. 



SIMULATION 

The previous analytic work tells us where to look if 
we wish to observe the terms in the oscillation formulae 
which are linear in sin 0x3. We investigate this further 
by utilizing the analysis from Ref. We proceed by 

fixing all of the oscillation parameters except 0x3 to their 
values that are given by an analysis in which 0x3 is set to 
zero. 

In Figs. 121 II we present the oscillation probabilities 
Vee, V e ^i, and V^p as a function of L/E respectively. In 
all cases we assume a Gaussian spread in energy of twenty 
percent. In all cases curves are presented for 0x3 = 
and 0x3 = ±0.2. In Fig. [2] for V ee there are only two 
curves as V e e is a function of sin 2 13 . Thus the curves 
for 13 = ±0.2 are identical. The curve is presented for 
completeness and to note that there is a measurable de- 
pendence of Vee on the magnitude of 0x3 near the peak at 
L/E = 3.2xl0 4 m/MeV, near the location of KamLAND 
as has been indicated in Ref. |l3j |. 

In Fig. |3 we verify two facts as derived analytically in 
the previous section. First, there is a significant linear 
dependence of V e)1 on sin 0x3 . Secondly, the dependence 
is maximal at tp\2 — tt/2 and 3 7r/2. The optimal value 
of L/E to measure 0i3, including its sign, is thus L/E = 
27r/A2x- The linear term in 0x3 clearly dominates near 
the maximum of the oscillation. However, the constant 
term in Eq. Q proportional to sin 2 0x3 becomes relevant 
whenever the probability is minimal. 
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FIG. 3: V ef i versus L/E. The solid curve corresponds to the 
mixing parameters for a fit with #13 = 0. The dashed (dot- 
dashed) curve corresponds to #13 = 0.2 (—0.2). 

In Fig.0]we see that is a somewhat better quantity 
to measure than "Pen, as here the linear term is even 
more dominant. Again, the linear dependence on #13 is 
maximal at y> 12 = ir/2 and 3n/2. Fig. ^ demonstrates 
the inadequacy of parameterizing oscillation parameters 
as a function of sin 2 # 13 . Figs. El and ^demonstrate this 
more dramatically as the linear term in #13 is found to 
dominate the quadratic term for V efl and in this 
region of L/E. 

The ideal measurement would be to determine T 3 ^ 
over a range for L/E from approximately 1.6 x 10 3 
m/MeV to 1.6 x 10 4 m/MeV. The lower end of this range 
provides an overall calibration point where the affect of 
non-zero #13 is small while the upper end is the point 
where the effect is largest. To avoid the further compli- 
cations of the Earth's MSW effect, the energy should be 
less than about 100 MeV, the energy of the MSW reso- 
nance in the mantle of the Earth. For an upper limit on 
the energy of 50 MeV, this sets an ideal value L — 80 km 
and a range in energy for the muon neutrinos of 5 to 50 
MeV. For a smaller value of L, the range of the neutrino 
energies would have to be proportionally smaller. More 
likely one's neutrino source would produce neutrinos with 
energies of at least a few hundred MeV's. Here, matter 
effects become significant though the qualitative features 
of the curves in Figs. JSJ and (0J can be salvaged. The 
most notable changes for oscillations through the mantle 
in this scenario are the increased frequency of oscillations, 
necessitating a suitable change in the baseline, and the 
adjustment of the matter mixing angle 0^. However, the 
difference between the matter mixing angle 6$ and the 
vacuum value of the mixing angle can be minimized by 
choosing a neutrino energy twice that of the resonance 
energy. For such a situation, the qualitative features of 
the oscillations presented above remain the same. 

Finally, in Fig. [SI we present the values of V ee , V eil: and 
Vpp as a function of #13. The other parameters are fixed 
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FIG. 4: The same as Fig. except the muon survival proba- 
bility 7-V M is presented. 

at their values for a fit with #13 = 0. We use a value 
of L/E = 1.6 x 10 4 m/MeV chosen so as to maximize 
the relative importance of the linear terms in #13. Again, 
Vee is quadratic in #13. However, V ef i and are nearly 
linear in $13 over this rather large range of —0.4 < #13 < 
0.4. The near linearity reinforces our observation that 
the measurement of V efJl and "P MAt at L/E — 271-/A21 is a 
way of determining #13, including its sign. 

The phenomenology of neutrino oscillations, in the ab- 
sence of CP violation and various exotica such as a fourth 
sterile neutrino, has been performed in the context of 
determining three mixing angles and two mass-squared 
differences. Historically, the results were presented in 
terms on sin 2 #13, which yields an upper limit. When 
the bounds on the mixing angles were explicitly quoted, 
they were all stated to be bounded by ir/2. In Ref. |2lj . 
it was shown that, in the physical case which necessarily 
includes the MSW matter effect, a second branch corre- 
sponding to 5 — 7r must also be included. In Ref. we 
extended the derivation, again in the physical case which 
necessarily includes the MSW matter effect, to show that 
only the 5 = branch is required if the mixing angle #13 
is allowed to vary from — ir/2 to +ir/2. There are two ad- 
vantages to this convention. First, the allowed region for 
6*i3 consists of a single region that extends on either side 
of 0, rather than two disjoint regions, one for 5 = and 
one for 6 = ir. Secondly, in the presence of CP violation, 
the CP phase is bounded by and ir. Thus a measure- 
ment that is sensitive to cos S could uniquely determine 
the quadrant in which S lies. 



CONCLUSION 

The natural question that arises from this formal work 
is whether you can measure #13 including its sign. In 
Ref. 01 we demonstrated that for a model analysis of 
the world's data, the % 2 space was asymmetric in #13, 
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FIG. 5: The oscillation probabilities versus the mixing angle 
9i3 with the other parameters fixed at their value from a fit 
with #i3 = 0. The solid curve is V ee , the dashed curve Ve^, 
the dot-dashed curve We take L/E = 1.6 X 10 4 m/MeV 
where the effect of the linear terms in O13 are maximized. 

thus demonstrating that such a measurement might be 
possible. A more recent analysis 13) which includes re- 
cent data verified this. Here we address the question di- 
rectly and demonstrate that the answer to our question 
is yes - the mixing angle #13 does appear linearly in the 
oscillation probabilities at a level where it dominates the 
quadratic term for the correctly chosen oscillation prob- 
abilities when measured at the correct value of L/E. 

We find that this is true for measurements of V eil and 
Vfifi at L/E = 271-/A21. This fortuitously also corre- 
sponds to a value of L/E where the contribution from 
the CP violating effects is minimal. A consequence of this 
work is that it reinforces our earlier thesis that parame- 
terizing neutrino oscillation probabilities as a function of 
sin 2 0i3 is inadequate. 

There remains additional work to be done. We find 
that present data indicate a correlation between the al- 
lowed value of #13 and #23 when the linear terms are 
included in the analysis. We are investigating the impli- 
cations for #13 that will result from measurements of #23- 
We have assumed no CP violation in our model analysis. 
Since the CP phase and #13 are interrelated, this needs 
further clarification. For cleanliness of interpretation, we 
propose do do experiments in a region of E where the 
Earth MSW effect is small. We are examining in a quan- 
titative way how the Earth MSW effect might modify an 
analysis should this be necessary. We, like others, have 
excluded the LSND Q experiments. If MiniBoone veri- 
fies the LSND results, then a whole new physics will be 
needed to reach a consistent understanding of neutrino 
oscillations. 
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